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Abstract. 

The constraint of incomprcssibility is often used to simplify the magnetohydrody- 
namic (MHD) description of linearized plasma dynamics because it does not affect 
the ideal MHD marginal stability point. In this paper two methods for introducing 
incomprcssibility are compared in a cylindrical plasma model: In the first method, 
the limit 7 — * 00 is taken, where 7 is the ratio of specific heats; in the second, an 
anisotropic mass tensor p is used, with the component parallel to the magnetic field 
taken to vanish, p\\ — > 0. Use of resistive MHD reveals the nature of these two lim- 
its because the Alfven and slow magnetosonic continua of ideal MHD are converted 
to point spectra and moved into the complex plane. Both limits profoundly change 
the slow-magnetosonic spectrum, but only the second limit faithfully reproduces the 
resistive Alfven spectrum and its wavemodes. In ideal MHD, the slow magnetosonic 
continuum degenerates to the Alfven continuum in the first method, while it is moved 
to infinity by the second. The degeneracy in the first is broken by finite resistivity. 
For numerical and semi- analytical study of these models, we choose plasma equilibria 
which cast light on puzzling aspects of results found in earlier literature. 
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1. Introduction 

We devote our attention to the ideal and resistive MHD models, which despite their 
dramatic simplification of plasma behaviour, are crucial to the design and operation of 
controlled fusion devices and are at the core of many astrophysical plasma models. 
Simplified models such as these have utility if they can be used to make testable 
predictions or if they yield insight into the internal processes of a system. We test 
two variants of incompressible resistive MHD with these criteria in mind. The starting 
point for the analysis of most plasma models is an understanding of the wavemodes that 
arise: this provides information about the linear response and stability of the system 
and provides a basis for much nonlinear analysis. We focus on the linear behaviour 
of the plasma in this paper. For the resistive MHD model, which includes dissipation, 
the wavemodes are non-normal: a full picture of linear plasma behaviour requires an 
analysis of the transient behaviour of the system, as well as the eigenvalue analysis 
which predicts asymptotic behaviours over long time scales. These are closely related 
via pseudospectral methods [T] . In this paper we restrict attention to eigenvalue analysis. 
For many plasmas of physical interest it is true that the resistive term is small: typically 
this is quantified by a large magnetic Reynolds number. It might have been expected 
that for small enough resistivity, the resistive MHD model could simply be treated 
as a perturbation of the ideal MHD model. However, the change induced is actually a 
singular perturbation, which introduces higher spatial derivatives. One of the interesting 
effects of this property is that eigenfrequencies in the ideal model are not necessarily 
approached by the eigenfrequencies of any resistive modes, even for vanishingly small 
resistivity. 

Many papers have been published on the stable resistive MHD spectrum and several of 
the early papers (j2J- 0) focused on cylindrical models. These papers have established 
certain generic features of the resistive spectrum. The resistive spectrum is discrete, 
unlike the ideal MHD spectrum which has continua: on some intervals, every frequency 
corresponds to a generalised wavemode. In the resistive spectrum, a large number of 
fully complex eigenfrequencies can be found, and in general these lie along loci, or lines, 
on the complex plane. Generally as the resistivity is decreased to zero these lines become 
densely populated with eigenvalues. 

Resistive MHD is a simple closure of the full kinetic equations, and as a result the plasma 
dynamics parallel to the magnetic field lines are often quite poorly represented [Zj. For 
Alvenic modes, which do not strongly compress the plasma, these parallel dynamics 
are generally unimportant. However, for the slow and fast magnetoacoustic waves, the 
parallel dynamics and the effects of compressibility are important; these waves are not 
necessarily well modelled by resistive MHD. 

It is possible to find the compressible resistive MHD spectrum numerically (as in £>]) and 
ignore the slow and fast magnetoacoustic waves that are present. On the other hand, 
there are approaches which promise to isolate the Alfvenic portion of the spectrum and 
simplify the analysis. We present two of these incompressible approximations, in which 
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the predicted motions of the plasma satisfy V-v = (at least approximately). One 
approach is to artificially set the ratio of specific heats 7 to infinity (as in [1] and 
In the other, an anisotropic mass tensor p is used, with the component parallel to the 
magnetic field taken to vanish, p\\/p± — > 0. With this density tensor, ideal eigenmodes 
are incompressible, but to ensure exact incompressibility for resistive eigenmodes 7 
must again be set to 00. We can view these models as the extreme cases of a generalised 
resistive MHD model with two parameters, 7 and p\\/p±- The two extreme cases are 
not equivalent, and the resulting spectra are qualitatively different. We investigate 
these two methods, and compare them with the compressible resistive MHD model. We 
specialise to equilibria with zero background flow. Note that 7 — > 00 may be physically 
appropriate for particular conductive fluids and plasmas with (3 » 1. 
First, we examine the plane waves of the homogeneous incompressible MHD model. 
Then we evaluate spectra in a simple cylindrical equilibrium for varying values of 7, 
and with and without an artificial anisotropic density. This illustrates the transition 
between the compressible and incompressible cases. We then discuss the spectra of more 
general plasma configurations. A WKB analysis of a generic incompressible model is 
then undertaken in order to understand the features of these spectra and to verify the 
numerics. We begin by solving the dispersion relation. Then the singular features of this 
function are explored by reducing it to a simpler form. To complete the groundwork for 
semi-analytic calculations, the behaviour of the wave equation near these singular points 
is examined. Finally, we use our WKB analysis to find the spectrum of an example case. 



2. Wavemodes in incompressible MHD limits 

The first step in the analysis of these incompressible limits is a determination of the 
wavemodes in a simple homogeneous plasma. To this end we follow [Zj and derive wave 
frequencies. We begin by considering a wave with wavevector at some angle to the 
magnetic field B = Bqz, so k = knz + k±x, travelling in a plasma with sound speed 
K — (7Po/Po) 1//2 an d Alfven speed V a = (B 2 /fi p Y^ 2 . We recover the Alfven spectrum: 



.2 _ ;„2ta2 
A 



HVa\ (1) 



and also two other solutions to the plasma equations: 
I p||fc 



where 



In low-/? compressible plasmas, oj + corresponds to the fast magnetoacoustic wave, and 
co>_ to the slow magnetoacoustic wave. 

In the limit 7 — > 00 (with p\\/p± = 1) we find uj\ — > 00 and uj 2 _ — > k 2 V a 2 , so that the slow- 
mode is now degenerate with the Alfven mode. In more general plasma configurations, 
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the slow and the Alfven wavemodes still occur at very similar frequencies, and therefore 
can be strongly mixed. We show this does occur, so that generic spectra determined 
are composed of an unphysical combination of these types of wavemodes. In the limit 
P\\/p± — » we again have uj + — > oo, but oj 2 _ — > k 2 V 2 + kj_V 2 , which is slightly larger 
than the fast magnetoacoustic frequency. In this case we have suppressed the slow 
magnetoacoustic waves. 

If we set p\\/p± — > 0, we can show from the linearised equations that for general resistive 
MHD wavemodes p\\ — > implies: 

Bo .V,V.v) =-^M. (4) 

where B and Pq are the equilibrium field and pressure, r\ is the resistivity and j and 
v are the perturbed current and velocity. So for the ideal case (77 = 0) we have that 
V-v is a constant on all irrational surfaces, and, by continuity, for finite toroidal or 
poloidal mode number, we must have V-v = 0. In the resistive case, we have small 77, 
but possibly large d/dr so that resistive modes are not strictly incompressible. However, 
if we also require 7 — ► 00 then the resistive modes are strictly incompressible. 



3. Numerical results of varying incompressibility 

In order to show the effect of incompressibility on the resistive MHD spectrum, we 
solved the compressible, resistive MHD equations numerically. We use a code based on 
the description in 

We examine a cylindrical, zero-shear model described in P , with 

(3 « 4%. The incompressibility is explored by varying 7 in the range 1 - 1000. The 
incompressible limits correspond to 7 — > 00, but in this case 7 ~ 1000 is high enough 
to demonstrate the limit. We define the magnetic Reynolds number S = tr/t^ where 
T4 and tr are the Alfven and resistive timescales. For a cylinder of radius r p we have 
t~a = rp(pop) 1 ^ 2 1 B z , and tr = r^Uo/^o- The magnetic field perturbations are of the form 
b = exp (imO — inz/r p — iut) b (r), with k = nr p /R, by analogy with the toroidal case 
R can be interpreted in the sense that 2irR is the length of the plasma column and n 
is the 'toroidal' mode number. We have j3 m 4%, which allows the slow-mode spectrum 
to be shown on the same scale as the Alfven spectrum in the compressible case. The 
resulting spectra are shown in figure d For this case m — 1, k — 0.35, nq = 1.2 and 
S — 1 x 10 4 . Note that figure |l(e)| corresponds to the limit 7 — > 00 and figure |l(f)| 
corresponds to the limit p\\/p± — > 0. 

In these cases the ideal Alfven continuum degenerates to a point, at u>a = 0.057, but the 
ideal slow continuum is finite in extent because of pressure and field strength variation 
across the plasma. The slow continuum extends to the origin because the pressure is 
taken to be zero at the plasma boundary. Note the fork structure seen for the slow modes 
near the origin of figure |l(a)| This fork structure is lost as 7 is increased [figures |l(b)f - 
(e)]. Finally, as 7 — >• 00, most of the mode frequencies are in the vicinity of a semicircle 
of radius uja on the complex plane. From the figure, we see that there are many more 
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Figure 1. The compressible resistive spectrum of a constant current model [J3 = 4%) 
for various values of 7. The ideal slow-mode continuum is represented by a grey line 
on the real axis. 



modes near uja in the 7 — > 00 model, than in the more physical compressible model. It 
has been shown in |Hj that for this incompressible case wavemodes are eigenfunctions of 
helicity and none of the modes correspond directly to physical compressible wavemodes. 
In figure |l(e)| the two loci of eigenvalues correspond to wavemodes of opposite helicity. 
For the p\\/p± — > model, we find a spectrum [figure |l(f)| very similar to the 
compressible spectrum in figure [T(a)j but with the notable absence of the slow-mode fork. 
The position of individual Alfvenic eigenvalues is in fact well preserved in this model. 
The only noticable deviation is the eigenmode near the real axis, at Re(oS) ~ 0.035, 
which has a frequency shift of magnitude ~ 0.004 as a result of setting p\\/p± — > 0. Since 
this Alfven eigenmode is fairly close in frequency to the slow modes, it is not surprising 
that it is the one most strongly modified by an assumption of incompressibility. 

4. Generic spectra in resistive MHD 

For general plasma configurations with shear, the resistive Alfven spectrum is usually 
found to form a fork (e.g. figure 121 or those in [2] - The rather different shape of the 
spectral loci in [fig ures[gg}(f)] is a consequence of the equilibrium having an Alfven 
spectrum which degenerates to a point. 

The fork structure in the resistive MHD spectrum has been qualitatively explained in 
terms of WKB analysis by examining turning points within the plasma, see [Sj and 
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Figure 2. A typical eigenvalue spectrum for a compressible (7 = |) resistive MHD 
case showing the complex frequencies of damped and growing normal modes. In this 
low-pressure case, slow magnetosonic modes have eigenfrequencies very close to the 
origin and are not shown. The plasma model parameters are similar to the constant 
current case, but with a small shear: m = 1, K = 0.35, nq(r) = 1.2 X (1 — O.lr) and 
S = 3 x 10 3 . The ideal Alfven continuum is represented by a grey line on the real axis. 




The fork has three lines joining at a point below the ideal MHD continuum. Two 
lines run between the intersection point and either end of the Alfven continuum. The 
third line runs around approximately in a quarter circle to touch the imaginary axis. 
In a simple model with toroidal current density constant across the plasma, there is 
an analytical solution for the 7 — > 00 resistive MHD spectrum ||. We show that a 
perturbed variant of this constant current model, in which a slight shear is given to the 
magnetic field, is still amenable to the manipulations performed in [0]. By introducing 
shear we produce a model which has a finite width Alfven continuum, in which we might 
hope to recover the generic fork structure found in compressible results. We therefore 
solved this model using WKB analysis to explain the qualitatively different spectrum. 
In the remainder of this paper we set p± = pn = p. 

5. WKB analysis of a small shear equilibrium in the limit 7 — ► 00 

The model case is derived from [0], which considers a cylindrical plasma with a constant 
axial field and no shear. This model has been studied earlier in JU|. The equations 
used for this analysis are those of linearised, resistive, incompressible MHD, with 7 — > 00: 

P/io^(Vxv) = Vx(B-Vb + b-VB), (5) 

and magnetic field given by Ampere's law 

? = Vx(vxB)-Vx(^-Vxb). (6) 

at po 

The curl of the equation of motion is taken in order to suppress the perturbed pressure. 
Also, we specialise to an equilibrium state with no plasma velocity. 
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The idea is to introduce the shear as a small quantity, of the same order as the inverse 
wave number. The analysis is then the same as the shear-free case, up to two orders in 
the inverse wavenumber. The radial dependence is included in the dispersion relation 
in the radially varying quantities: B z (r), B p (r) and q(r). We first take the large 
wavenumber limit by ordering V ~ 0(l/e). For significantly dissipative modes, maximal 
balance of Ampere's law (JHJ) occurs for e ~ 0(S~ l l 2 ). In a typical physical situation we 
might have S > 1000 and thus e < 0.03 is a good expansion parameter. 
The magnetic field is expressed as B = zB z (r) + r6B p (r) with d(log[B p (r)}) / dr and 
d(log[B x (r)])/dr both of O(e), in order to satisfy the requirement of small shear. We 
again look at perturbations of the form b = exp (imd — inzjr v — iut) b (r) . For 
convenience we set b as 0(1) and this then implies v to be of 0(1) to complete the 
ordering. By using the relations V • b = V • v = 0, equations (J5J) and (0) can be reduced 
to: 

p^o— (V x v) = ^-[m - ng(r)]V x b + ^fi-nh + O(e) (7) 

and 

- — b = ^M[ m _ nq (r)}v - ^-(V x V x b) + O(e). (8) 
t~a r p p 

The safety factor q(r) is given by r p B z (r)/ RB p (r) and the non-dimensional resistivity 
r] = S- 1 . 

In this form, the only differential operator is the curl operator. This motivates us to 
look for solutions which are eigenfunctions of this operator, suggesting the ansatz 



r p V x v = av, r p V x b = ab, (9) 

which solves equations ((7|) and (jHJ) provided 

/ \ 2[m-nq{r)]n 

a(r,u) = - )2 . (10) 

B zy -^ 2 iuj [iu> — S^a^r, uo) 2 ] + [m — nq{r)] 2 
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By taking the curl of equation JOJ) we get 

a 2 

V x V x v = -V 2 v = — v, (11) 

r p 

since the velocity is divergence-free. This implies a relation for the z component of v 

1 d d f a2 + 2 m 2 \ 

r dr dr z \r 2 r 2 J z 

This is amenable to standard WKB analysis if a is large, and in this WKB limit equation 
(JT2J) is equivalent to: 

|U = -Q(r)v„ (13) 

with Q(r) = a 2 /r 2 . This will break down near the origin (r = 0) where we will use a 
Bessel function matching. Equation (jlHj) is solved approximately by: 

v, ~ a out Q-^ + a in Q-2e-^ (14) 

where the amplitudes a ou ^ and slowly varying functions, and 

/r 
Q^(r')dr' (15) 

Thus a/r p is the radial wavenumber and equation fTUj) provides the dispersion relation. 
6. Characterising the Stokes points 

To find the WKB solutions, it is first necessary to examine the structure of the 
dispersion relation in the plasma region. In particular, singularities and zeros and 
the associated branch structure of the dispersion relation must be examined. Branch 
points of the dispersion relation are known as Stokes points. The dispersion relation 
(fTUJl can be written as a cubic equation in a, with the coefficients as functions of 
q(r) = r p B z (r)/RB p (r), i.e. 

a 3 U^^] +a([m-nq{r)} 2 - (^) u 2 I = 2[m - nq(r)] K . (1(3) 



We would like to discover the singularity structure of our dispersion relation. Solving 
equation ()16|) for a leads to very ungainly equations and proves not to be enlightening, 
so we look for a simpler relation which will be topologically equivalent. Let us consider 
the case where there is no magnetic surface resonant with the perturbation. In this case 
we have [nq(r) — m] ^ within the plasma, and assuming also q(r) ^ then we can 
divide through the equation by 2[m — nq{r)\K and introduce a new variable a so that 

a 3 + x(r)a = l, (17) 
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a = a- 



and 



x(r) = 



{2[nq(r) — m]KS}3 
{q{r)R/r p )i 

{[nq(r) — m] 2 — [qR/r p ] 2 tu 2 } 



{2[nq(r) — m}nS(qR/r p )} 3 
The solution of equation (|17j) for a is 



;is) 



(19) 



23 x 



(27 + V729 + 108 x 3 ) ; 
where ^ is one of the cube roots of —1: 



(27 + V729 + 108 x 3 ) 1 



323 



(20) 



l + zv^ l-iV3 



(21) 



2 ' 2 

We consider i as a new radial variable. The function a is represented graphically by 
the Polya plot in figure EJ 





Figure 4. One branch of the multivalued function a(x) (£ = —1) shown on the 
complex plane as a Polya plot. Branch cuts are indicated as thick lines. The vector 
(Re [a], — Im[a]) is displayed on a grid. 

By inspection of the form of equation ()20j) we have candidates for branch points at the 
three roots of a; 3 = —27/4. However, not all of these candidate branch points are true 
branch points, as suggested by figure 0] This can be seen in the £ = —1 case, where we 
have, from equation fTTj) 
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which can be considered as a local inverse of equation (J2UJ) . 

Let us consider the neighbourhood of x = — ^/27/4 (using the principal value, so x is 
a negative real). We might expect a branch point here from the structure of equation 
(I20J) . At x = —^27/4 we have a = 2 2 / 3 . However, there is a neighbourhood around 
a = 2 2 / 3 where equation (|2*2*|) is analytic and has a non-zero derivative, and therefore 
the function has an analytic inverse around x = — y/27/4. There obviously cannot be 
a branch point in an analytic region. The other two candidate points are true branch 
points. It similarly follows that each of the other cases of equation (|2()jl have only 
two branch points each. Note that around a Stokes point at some position xq, we do 
not have a oc \Jx — Xq, as is typical for many WKB analyses [11J. Instead, we have 
a ~ C + Dyjx — Xq. 



7. Phase matching: a solution near the singularities 

In order in proceed with WKB analysis, we need to determine the behaviour of 
solutions near the Stokes points, the branch points of the dispersion relation. In the 
neighbourhood of the branch point, we approximate the dispersion relation by: 

Q(x)~l+Ax^. (23) 

This is unlike the more usual situation in WKB analysis where Q(x) ~ x around the 
Stokes points. The simplest treatment of the phase matching follows from considering 
A <C 1 in which case the A = case can be used as a zeroth order solution in a region 
around the Stokes point. Note that for A = 0, the dispersion relation is independent of 
x and there is no reflection of the wave. As we will see, as A — > the reflectivity goes 
to zero. The transmitted part of the wave will be decaying for finite A, so that we have 
partial absorption of the travelling wave. 
Our wave equation is 

d 2 y , , i, j i 

1 ± = -(1 + Ax^y t (24) 

with an A = solution 

Vo = e ±ix , (25) 

which motivates the substitution 

The other choice of sign leads to a second solution to the equation, which is growing for 
x — > — oo. Substitution of equation (}26|) into equation (Ell) leads to 



, r- du ( du\ 2 d 2 u , 
AV *- 2i Tx + [Tx) + ^ = °- (27) 
We are looking for small departures from the A = solutions and in this case we can 
choose du/dx 1 so that to first order 

, ._ du d 2 u , , 

A ^~ 2l Tx + dx^ = ^ (28) 
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from which we can find u'ix) 



u'(x) = e 2ix C- l -A 



4y£ + e (2 *^ /4)i V^erf (e"^ /4 v^x)]. (29) 



The coefficient of integration, C, must now be chosen such that we can match the 
solution on the left-hand side of the origin to the evanescent WKB solution. We have 
required u'{x) C 1, so an oscillatory u(x) can be modelled as ee 2lx with e< 1 (plus a 
constant which can be safely ignored) in which case: 

These correspond to the WKB solutions, which are approximately of the form c\e~ lx + 
C2e lx near the origin. We require that the WKB solution matched on the left-hand 
side have C2 = because the corresponding term grows exponentially for large negative 
x. We therefore must have e — > as x — > — oo. Using the asymptotic expansion of 
erfc = 1 — erf, as given by equation 7.1.23 of we find the x — > — oo limit of equation 
([29)1 . allowing us to express this matching condition as: 

C = ^(l + i)V^. (31) 
Then we have a solution for y which is asymptotically of the form: 

y = P(x) (e~* x + e ^v^(l-*n (32) 



for x ^> 1, with P(x) a slowly varying function. The phase matching condition is given 
by finding the nodes of these waves, which fixes the WKB phase at x = 0: 

; ~i i Al+i)Ay/7T 

0o = -j log( ). (33) 

8. Finding wavemodes 

Global modes are found in the usual way: we look for paths C in the complex plane 
joining the axis and boundary where J B a(x)dx is real for any sub-path B of C. These 
paths will be WKB solutions if the integral f c a(x)dx = Jj - a(x)dx which can be 
guaranteed if there are no singularities of our differential equation coefficients in the 
region. In particular, this requires that the circular path C — [0, 1] does not enclose any 
Stokes points. The quantisation condition is supplied by requiring the correct behaviour 
at boundaries. At the origin the WKB wavemode must be matched to a Bessel function, 
and this gives the condition 4>\ x= q = 7r(l/4+m/2). At the outer boundary of the plasma, 
we require v r — (fixed plasma boundary), which leads to <f>\x=i = tt/2. 
Localised modes proceed from the axis or outer boundary of the plasma and propagate 
along ray trajectories (which will in this case be anti-Stokes lines) to a Stokes point. 
They are then evanescent past this point, so it must be possible to draw a path 
connecting the relevant Stokes point to the other boundary without crossing a Stokes 
line. For Stokes points of the form Q(x) = a + bx^, which are present in this analysis, 
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we have a complex phase matching criterion. The phase integral between the Stokes 
point and the boundary is then required to be complex for matching to occur. This 
means that we cannot follow anti-Stokes lines, along which the phase is real, exactly to 
join the boundary and the Stokes point. The complex portion of the phase leads to a 
correction to the path, which must be taken into account. 



9. Application of the WKB method to the small shear incompressible case 

For explicit studies, we use a small shear test case: 



BJ r ) = = 2?4 , (34) 

pK J 1-rS 1-O.lr' V ; 

20 

B z {r) = 1, R = — , n = 10, m = 1, r p = 1, q{r) = 0.12(1 - O.lr) (35) 

The WKB trajectories in the complex r plane were determined numerically, and several 
loci found by finding appropriate paths in the plane, as in [21 • The loci can be 
characterised by the branch of the dispersion relation which they lie on, and whether 
the corresponding wave modes are fully global modes or have a turning point inside the 
plasma. 

Although there are three branches of the dispersion relation, on one branch it is never 
possible in practice to form global modes: the rays inevitably escape towards complex 
infinity. The other two branches then produce the two forks. 

The eigenvalues are displayed in figure El together with the numerical result from a code 
based on jH]. The spectrum is qualitatively similar to a fork structure, but also shares 
the features of the original simple model. Note that the double loci (running parallel to 
each other in an arc) are still present in this model. 
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Figure 5. The resistive spectrum from numerical solution of the MHD equations 
(circles) compared to the WKB result (plus signs). The ideal Alfven continuum is 
represented by a grey bar on the real axis. Plasma parameters are given by equation 
(ESJ and S = 3 x 10 4 . 
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The nature of the difference between the two branches of the double locus can be seen 
in equation ©, and the form of the dispersion function for large enough a. Here we 
have two solutions for a(u) such that a± ~ — ao, and the two WKB solutions consist of 
waves of opposite helicity. Finite pressure gradients in this equilibrium result in waves 
of opposite helicity having slightly different frequencies. 



10. Effects of the 7 — > 00 approximation 

The reason why we see a pair of loci in l(e)| rather than the single locus usually depicted 



for compressible spectra (e.g. figure 1(a) ) is that in this incompressible model (the limit 
7 — ► 00) there are two classes of wavemodes present which can be excited at the Alfven 
frequency. In a uniform field these wavemodes are degenerate: they oscillate at the same 
frequency. However the two frequencies are split when the plasma contains currents 
perpendicular to the magnetic field (i.e. in non- force- free plasmas). In the compressible 
model at low /3, these two degrees of freedom correspond to the slow (magnetosonic) 
and Alfven wavemodes and the ratio between slow frequencies and Alfven frequencies 
is of order /5 1//2 . 

Force-free models are important special cases, in which pairs of loci of eigenvalues 
coincide. The 7 — > 00 approximation will still result in unphysical eigenmodes. We 
note the paper of Ryu and Grimm jl], which uses this incompressibility assumption 
to analyse a case with finite pressure gradients where we should see a double locus 
structure. We nevertheless see a simple fork structure. We believe that the splitting 
effect is rather small in this case, so that what looks like one fuzzy locus is in fact a 
double locus. 



11. Conclusions 



In plasma physics an assumption of incompressibility is often justified because the 
parallel dynamics of the plasma and the fluid compression across the field are much 
less important than the forces due to the magnetic field. For example, incompressibility 
does not generally affect ideal MHD marginal stability (but this does not extend to 
resistive MHD [13]). 

Two incompressible resistive MHD models were compared with the physical 
compressible model by analysis of their spectra. For the first model, where the ratio 
of specific heat is taken to infinity, we expect from local analysis to find two types of 
wavemodes present at the Alfven frequency. In the second model where we again set 
7 — > 00, the parallel plasma inertia is set to zero, and we expect only one Alfvenic 
mode to be present in the spectrum, corresponding to the physical case. Numerical 
computation of the spectra of a magnetic shear-free plasma confirms that the second 
model reproduces most of the eigenmodes associated with the Alfvenic model correctly. 
The first model has twice as many modes present at the Alfven timescale. 
It is noted that in general most of the modes resolved do not correspond to Alfven 
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modes and have no physical significance. The shape of an incompressible spectrum for 
a more general model, with shear present, was determined numerically and by WKB 
analysis. The unusual nature of the local dispersion relation leads to a complex structure 
of loci. The resulting spectrum included the 'double locus' of the zero shear model and 
also demonstrated the fork structure that is seen generically for stable resistive MHD 
spectra. 

There are many qualitative features of the resistive Alfven spectrum that can be 
reproduced by simply setting 7 — > 00. Unfortunately, physical wavemodes and 
frequencies are not well modelled in this approximation. The stable part of the ideal 
Alfven spectrum is irreparably mixed with spurious modes in this limit. However, 
by using an anisotropic mass density tensor, an incompressibility constraint can be 
introduced while preserving the Alfven modes. 
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